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ABSTRACT 

The self-binding energy and stability of a galaxy in MOND-based gravity are curiously decreasing 
functions of its center of mass acceleration (of the order of 10 -12 - 10 _10 m/s 2 ) towards neighbouring 
mass concentrations. A tentative indication of this breaking of the Strong Equivalence Principle in 
field galaxies is the RAVE-observed escape speed in the Milky Way. Another consequence is that 
satellites of field galaxies will move on nearly Keplerian orbits at large radii (100 - 500 kpc), with a 
declining speed below the asymptotically constant naive MOND prediction. But conseqences of an 
environment-sensitive gravity are even more severe in clusters, where member galaxies accelerate fast: 
no more Dark-Halo-like potential is present to support galaxies, meaning that extended axisymmetric 
disks of gas and stars are likely unstable. These predicted reappearance of asymptotic Keplerian 
velocity curves and disappearance of "stereotypic galaxies" in clusters are falsifiable with targeted 
surveys. 

Subject headings: gravitation - dark matter - galaxies: structure - galaxies: kinematics and dynamics 



1. INTRODUCTION 

The exact nature of dark matter is an outstanding puz- 
zle despite our ability to carry out increasingly realistic 
simulations. No astroparticle or gravitational theory so 
far can account for its various effects on both galactic 
and large scales satisfactorily (e.g., Zhao 2007). 

In relatively isolated (field) galaxies, observations of 
a tight correlation between the mass profiles of bary- 
onic matter and dark matter at all radii (McGaugh et. 
al. 2007; Famaey et. al. 2007a) are most often inter- 
preted as supporting the modified Newtonian dynamics 
(MOND; Milgrom 1983; Bekenstein & Milgrom 1984). 
Indeed, without resorting to galactic dark matter, this 
simple prescription reproduces (to amazing consistency) 
the kinematics of galaxies over five decades in mass (e.g. 
Sanders & McGaugh 2002; Bekenstein 2006; Milgrom & 
Sanders 2003; Sanders & Noordermeer 2007; Nipoti et 
al. 2007; Famaey, Bruneton & Zhao 2007b; Gentile et. 
al 2007ab). The recently devised relativistic counterparts 
of MOND (Bekenstein 2004; Zlosnik, Ferreira & Stark- 
Electronic address: hz4@st-andrews.ac.uk 



man 2007) also enabled to study the Cosmic Microwave 
Background (Skordis et al. 2006) as well as gravitational 
lensing of galaxies and galaxy clusters (e.g., Zhao et al. 
2006; Chen & Zhao 2006; Angus et al. 2007; Famaey et 
al. 2007c). However, it is also important to observation- 
ally and theoretically examine how internal properties of 
galaxies like their rotation curve and morphology should 
change with the environment. 

Indeed MOND-based theories generate different de- 
grees of dark matter-like effects depending on the abso- 
lute acceleration. Most galaxies are, like the Milky Way, 
in the field where they accelerate slowly with respect to 
the Cosmic Microwave Background, typically at a rate 
of 0.01a to 0.03o (Famaey et al. 2007b; Milgrom 2002; 
Angus & McGaugh 2007). But in X-ray clusters, galax- 
ies accelerate much faster, from 0.3ao to 3ao (Angus et. 
al. 2007; Pointecouteau, Arnaud & Pratt 2005). 

This external gravitational field has wider and more 
subtle implications for the internal system in MOND 
than in Newton-Einstein gravity, for the very reason that 
MOND breaks the Strong Equivalence Principle. In par- 
ticular, it is well known that MOND potentials are loga- 
rithmic for isolated distributions of finite mass, and con- 
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sequently infinitely deep, but that the internal potential 
becomes "polarised Keplerian" at large distances (Beken- 
stein & Milgrom 1984; Zhao & Tian 2005) when an ex- 
ternal field is applied. 

In this contribution, we numerically solve the MOND 
Poisson equation for systems embedded in several dif- 
ferent environments, ranging from the field to galaxy 
clusters, and show (i) that for the Milky way (embed- 
ded in a weak gravitational field), the local escape speed 
is numerically compatible with the observations as an- 
alytically predicted in Famaey et al. (2007b), (ii) that 
rotation curves of Milky Way-like galaxies would have a 
rapid Keplerian fall-off when residing close to the center 
of clusters, while this fall-off in field galaxies would hap- 
pen at 100 - 500 kpc, and (iii) that usual Low Surface 
Brightness disks should not exist in MONDian clusters. 

2. BINDING ENERGY OF AN ACCELERATING 
MILKY WAY 

Galaxies free-fall, but with slowly-changing systematic 
(center-of-mass) velocity v com (t). Their present non- 
zero systematic velocity is mainly the accumulation of 
the acceleration by the gravity from neighbouring galax- 
ies over a Hubble time. Consider, as a first approx- 
imation, that a galaxy is stationary in a non-inertial 
frame (in the Galilean sense) , which free-falls with a "uni- 
form" systematic acceleration v com = g cxt X = est due 
to an external linear potential, say, —g cyA X along the 
A-direction, where the over-dot means time-derivatives. 
Let Vi n t = (X,Y,Z) be the peculiar acceleration of a 
star-like test particle in the coordinates relative to the 
center of a non-evolving galaxy internal mass density 
p(X,Y,Z), then 



Vint = g 



= -V$ int (A,Y,Z), 



(1) 



where g is the absolute acceleration g satisfying the 
MOND Poisson's equation 



V.I/x^g] = 4mGp{X,Y,Z), 



(2) 



One can define an "effective potential" $; nt (A, Y, Z) 
(called "internal" potential) and an "effective energy" 
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E c s = ^f 1 + <E>int(A, Y, Z), where E e ff is conserved 
along the orbit of the test particle effectively moving 
in a force field — V ( I > i n t(A, Y, Z), which is curl-free and 
time-independent because the absolute gravity g is curl- 
free, center-of-mass acceleration v com is assumed a con- 
stant, and the galaxy density p{X, Y, Z) is assumed time- 
independent. 

Far away from the center of the free-falling system, 
we have |v int | < |v com |, hence p -> p, m = ^(v com /a ) = 
est, and the equation reads (Bekenstein & Milgrom 1984; 
Milgroml986; Zhao & Tian 2005; Zhao & Famaey 2006): 



d 2 



dX 2 



AirGp/pr, 



(3) 



where Y, Z denote the directions perpendicular to the ex- 
ternal field A-direction, and A = [dhip/dhix] x= ^ com y aQ 
is a dilation factor (note that 1 < 1 + A < 2). So at 
large radii where the external field dominates, and the 
equation is linearizable, the potential satisfies a mildly 
anisotropic Poisson equation, and the solution at large 



radii 1 goes to 



$% t (X,Y,Z) = 



GMw 



p m ^(l + A)(Y 2 + Z 2 )+X 2 + s 2 



( 4 ) 

where we included a softening radius s, comparable of the 
half-light radius of a galaxy. Hence the internal potential 
$i n t is finite, and approaches zero at large radii. 

The escape speed at any location r in the system can 
then be meaningfully defined by 



= E eS 



v 2 sc (X,Y,Z) 



$ int (X,Y,Z). 



(5) 



Such escape speed is a scalar independent of "the path to 
escape" because the "effective energy" E c g is conserved, 
and a particle with E c g equal zero (the maximum value 
of 'I'int) w ih reach infinite distance from the system, and 
never return, hence will be lost into the MOND potential 
of the background (from which it cannot escape) . How- 
ever, equal escape speed contours across a disk galaxy are 
generally not axisymmetric, meaning the escape speed on 
opposite symmetric locations of the Galaxy differ. 

Hereafter, we numerically solve Eq.[2]using the MOND 
Poisson solver developed by the Bologna group (Ciotti, 
Londrillo, & Nipoti 2006); the results based on spherical 
grids are also confirmed with the cartesian grid-based 
code of the Paris group (Tiret & Combes 2007) with 
very different spatial resolutions. We program in the 
mass density of the internal system, solving the MOND 
Poisson equation as if it were isolated, except for re- 
quiring a boundary condition on the total gravity as 
-g -► <? oxt X - V<S>™ t (X,Y,Z) on the last grid point 
(A, Y, Z). Note finally that in our models hereafter, we 
use the parametric /i-function p(x) = x/(l + x), which 
fits well the rotation curve of the Milky Way (Famaey 
& Binney 2005), as well as external galaxies (Famaey et. 
al. 2007a; Sanders & Noordermeer 2007). 

We use the Besangon Milky Way Model (Robin et. al 
2003) to simulate High Surface Brightness (HSB) galax- 
ies. This model is a realistic representation of the Galaxy, 
explaining currently available observations of different 
types (photometry, astrometry, spectroscopy) at differ- 
ent wavelengths. The stellar populations included in 
the model are : a thin disk made of seven isothermal 
layers each having a different age, between 0.1 and 10 
Gyr; a 11-Gyr-old thick disk with a modified exponen- 
tial density law, a spheroid with a power law density, 
slightly flattened, a prolate old bulge modeled by a tri- 
axial density law. We logically take the dark matter 
halo away from our simulations. We then apply the 
MOND Poisson solver using 512 x 64 x 128 grid points 
where the grid points in the radial direction are chosen 



as r, = 50.0 tan 



(i + 0.5) 



0.57T 



512+1 



kpc. 



As a first application, the RAVE solar neighbourhood 
escape speed 544+4g km/s is well-reproduced by our fully 
numerical model galaxy (Fig. 1 for a typical external 
field of O.Olao), as analytically anticipated in Famaey 
et al. (2007b). When the direction of the external grav- 
ity changes, the escape speed changes in a narrow range 

1 throughout the paper, by "large radii" we mean a distance 
large enough to neglect the internal field </; nt -C v C om = ffext but 
small enough that the external field g cx t can be treated as constant 
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Fig. 1. — Model of a Milky Way-like galaxy in weak (O.Olao) and 
strong (2ao) external fields. Escape speeds in the disk plane for 
various field directions (solid and dotted) are compared with the 
error bar for the local escape speed measured from the RAVE sur- 
vey(Smith et. al. 2007). The predicted circular speed curves (dot- 
dashed lines) are also compared with data (diamonds) (Caldwell & 
Ostriker 1981). Note that for an external field of O.Olao, the transi- 
tion radius between MOND and Keplerian regime is at ~ 500 kpc, 
while it is at ~ 150 kpc for an external field of 0.03ao. 
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Fig. 2.— Similar to Fig.l, but for models of a NGC 1560-like LSB 
galaxy. Circular speed curves for g ex t = (no escape), g ex t = 0.3arj 
(middle) and g cxt = 2ag (bottom) are compared with the observed 
rotation curve of NGC 1560. 



[545, 558] kms -1 in the solar neighborhood. We more- 
over point out that satellites of Milky Way-like galaxies 
will move on nearly Keplerian orbits at very large radii 
(> 100-500 kpc depending on the external field strength), 
with a speed much below the asymptotically flat naive 
MOND prediction. 

We also model NGC 1560, a benchmark Low Sur- 
face Brightness (LSB) disk galaxy (Broeils 1992). We 
use an exponential stellar disk of 1.97 x 1O 8 M and a 
multi-Gaussian gaseous component of 1.07 x 10 9 M Q to 
match the observed baryon distribution(Broeils 1992). A 
MOND Poisson solver (Ciotti, Londrillo, & Nipoti 2006) 
with 256 x 64 x 64 grid points is applied, the radial grid be- 

(i + 0.5)-°- 5n 



mg n 



10.0 tan 



'256+1 



kpc. Again it is found 

that an acceleration of O.Olao is compatible with the ro- 
tation curve of NGC 1560, which is an isolated galaxy. 



3. FAST-ACCELERATING GALAXIES IN 
CLUSTERS 

Now consider boosting the Milky Way's systematic ac- 
celeration suddenly to match the environment in a galaxy 
cluster. Fig. [1] shows for an external field of 2ao, the es- 
cape speed of stars is much reduced, falling Keplerian- 
like 300 y/5 kpc/r km s _1 outside 5 kpc, where half of the 
stars and gas of the Milky Way are located. All dwarf 
satellites of the Milky Way, and outer disk rotating with 
200km/s would then barely be kept from flying away. 

In fact, the instantaneous hypothetical circular speed 
must also be lowered by the sudden boost of acceleration, 
and the outer galaxy (> 5 kpc) should exhibit a Keple- 
rian falling rotation curve (Fig.l). Outer disk stars and 
gas should enter elliptical or parabolic orbits of the same 
angular momentum if allowed to respond to a suddenly 
reduced gravity, and precess with a preferred direction of 
instantaneous systematic acceleration which thickens the 
disk. In any case, observing asymptotically flat rotation 
curves for purely baryonic Milky Way-like galaxies resid- 
ing in such environments would falsify MOND. Galaxy 
of same luminosity would have lower velocities, consis- 
tent with the observed trend with cluster Tully-Fishcr 
relation(Sanders & McGaugh 2002). 

Now let us consider suddenly boosting the accelera- 
tion of our benchmark LSB galaxy to v com = 0.3ao or 
i>com — 2ao (typical of outer and inner parts of galaxy 
clusters. The real orbit of a member galaxy would pass 
both regions at apocenter and pericenter respectively). 
The circular speed (Fig. [2] lowest dot-dashed curve) of a 
fast-accelerating LSB is very much reduced. All previ- 
ous disrupting effects are even more severe on an LSB 
galaxy with the escape speed falling as low as 50 km/s. 
Outer stars with original circular speed 50 — 80 km/s 
would enter parabolic orbits, and inner stars move out- 
wards on severely elongated non-planar orbits. Actually, 
the dynamics resembles a purely Newtonian disk with- 
out a round stablizing dark halo, meaning that the galaxy 
would become extremely bar-unstable (Mihos, McGaugh 
& de Blok 1997). Such an LSB would lose its MOND 
support and would be subject to strong distortions, even 
before the traditional tidal effect becomes important. 

4. CONCLUSION AND DISCUSSION 

The external field effect is a generic prediction of 
the modified gravity theories where the modification is 
acceleration-based and violates the Strong Equivalence 
Principle, e.g., the relativistic versions of the AQUAL 
theory of Bekenstein & Milgrom (1984). The effect is 
helpful (i) to allow high velocity stars to escape from 
Milky Way-like field galaxies, and (ii) to decrease the 
orbital velocities of their satellite galaxies at very large 
radii (100 - 500 kpc), contrary to the naive expectation 
of MOND without external field effect, i.e. that rotation 
curves should be asymptotically flat. In this respect, the 
data of Klypin & Prada (2007) will be very useful in 
the future to verify/falsify/quantify this effect thanks to 
detailed numerical modelling. 

On the other hand the internal dynamical structure 
of a field galaxy would transform when entering a clus- 
ter. Classical relations of field galaxies, such as the 
Tully-Fisher relation, the galaxy luminosity functions, 
the Hubble type distribution etc. are expected to mod- 
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ify strongly in clusters. The effects are most destructive 
for classical LSB galaxies; curiously their field counter- 
parts have been a legendary success for MOND in terms 
of well-fitted rotation curves. 

We thus argue that it would be extremely valuable to 
analyse the kinematics of a sample of HSB galaxies and 
search for LSB galaxies in nearby clusters using deep 
HI surveys. The study of a sample of galaxies would be 
needed because of the uncertainty of the determination of 
the real distance (as opposed to projected distance) of a 
galaxy from the cluster center. An obvious difficulty will 
be that cluster galaxies are Hi-deficient (Solanes et.al. 
2001). An example of such an HI database is the VIVA 
survey (VLA Imaging of Virgo in Atomic Gas (Chung, 
van Gorkom et. al. 2007). We also predict that a future 
detection of any undistorted HSB late-type disk galaxy 
near the center of a galaxy cluster would be extremely 
surprising in the context of MOND. A null-detection of 
thin LSB disks is predicted in clusters because they most 
probably have been turned into gas-poor dwarf ellipti- 
cals if not fully disrupted. Low surface density gas in 
galaxies also suffers ram pressure stripping while mov- 
ing in gaseous clusters; gas is easily stripped in the re- 
duced MONDian internal gravity, further reducing avail- 
able mass for self-gravity. 

Surely similar effects occur in the context of cored Dark 
Halos. Some simulations show that LSB disks and dwarf 
irregulars get harassed (Moore 1999) and transformed 
into dwarf ellipticals or ultra-compact dwarf ellipticals 
(Evstigneeva et al. 2007; Cortese et al. 2007) in the 
densest part of the cluster coinciding with the region 
where the external field is the highest. To our knowl- 
edge the properties of cluster disk galaxies (such as their 
Tully-Fisher relation) have not been extensively simu- 
lated. The important assumption of existing simulations 
is a large core for the CDM halo of the LSB; the ha- 
rassment becomes much less effective if the cluster mem- 
ber LSB started with a dense CDM cusp density (Lucio 
Mayer 2007, private communications). 

2 For similar reasons open clusters in the solar neighbourhood 
are predicted internally unbound in MOND. 



As for the gas-poor non-rotating dwarf spheroidals 
(e.g., Sextans), they are expected to have a central CDM 
density of ~ O.IMqpc -3 , a factor of 100 denser than 
the galaxy cluster, hence might survive the tidal ha- 
rassment in CDM. If in MOND a spheroidal object of 
M = 2 x 10 5 M Q , and half mass radius of s ~ 500pc is 
suddenly introduced into a galaxy cluster, it would have 
a central binding energy of only ~ — (1 + A)" 2 / 3 /!" 1 ~ 
(5kms -1 ) 2 , much less than its initial internal random 
motion energy 3x ( 10k ™ s — L ; hence is quickly dispersed 
(perhaps anisotropically) . 2 In short, any discovery of a 
sample of classical LSB galaxies in clusters would favor 
cuspy CDM, and falsify MOND or cored Dark Halos. 

The tidal harassment effect exists in MOND as well 
(Zhao 2005); cluster galaxies suffer from tides in addi- 
tion to the unique destructive effect of the external field. 
An even more curious distortion to the MONDian LSB or 
HSB disk happens when the disk is mis-aligned by an an- 
gle 9 m with the instantaneous direction of the external 
field, which generally changes amplitude and direction 
along the orbit of an LSB on time scales of 0.2-lGyr. 
The elliptical potential of Eq.(4) creates a differential 
force with a component normal to the disk, hence a spe- 
cific torque — r x V$. This causes differential precession 
of the disk angular momentum vector with an angular 
speed proportional to /i" 1 yj GM/r 3 A sin(20 m ) ; an LSB 
disk is likely shredded by one precession and a HSB disk 
is thickened. The precession, asymmetric dilation and 
reduction of inner circular velocity curves (cf. Fig. 2) are 
confirmed by N-body simulations using the code of the 
Paris group (Tiret & Combes 2007) in MOND, but are 
generally forbidden by Newtonian laws in the Dark Mat- 
ter framework. 



We thank Luca Ciotti, Pasquale Londrillo, Carlo 
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ler for helpful comments on an earlier manuscript. 
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